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Abstract 

We have analyzed magnetized equilibrium states and showed a condition for appearance of the 
prolate and the toroidal magnetic field dominated stars by analytic approaches. Both observa¬ 
tions and numerical stability analysis support that the magnetized star would have the prolate 
and the large internal toroidal magnetic fields. In this context, many investigations concerning 
magnetized equilibrium states have tried to obtain the prolate and the toroidal dominant solu¬ 
tions, but many of them have failed to obtain such configurations. Since the Lorentz force is a 
cross product of current density and magnetic field, the prolate shaped configurations and the 
large toroidal magnetic fields in stars require a special relation between current density and the 
Lorentz force. We have analyzed simple analytical solutions and found that the prolate and the 
toroidal dominant configuration require non force-free toroidal current density that flows in the 
opposite direction with respect to the bulk current within the star. Such current density results 
in the Lorentz force which makes the stellar shape prolate. Satisfying this special relation be¬ 
tween the current density and the Lorentz force is a key for appearance of the prolate and the 
toroidal magnetic field dominated magnetized star. 

Key words: stars: magnetars — stars: magnetic fields - stars: rotation 


1 Introduction 

Anomalous X-ray Pulsars and Soft-Gamma-ray-Repeaters 
(SGR’s) are considered as special classes of neutron stars, i.e. 
magnetars (Thompson & Duncan 1995). According to obser¬ 
vations of rotational periods and their time derivatives, magni¬ 
tudes of global dipole magnetic fields of magnetars reach about 
1 q14-i5q Recently, however, SGR’s with weak dipole mag¬ 
netic fields have been found (Rea et al. 2010, 2012). Their 
observational characteristics are very similar to those of ordi¬ 
nary SGR’s but their global dipole magnetic fields are much 
weaker than those of ordinary magnetars. It might be ex¬ 
plained by a possibility that such SGR’s with small magnetic 


fields hide large toroidal magnetic fields under their surfaces 
and drive their activities by their internal toroidal magnetic 
energy (Rea et al. 2010). Recent X-ray observation of mag- 
netar 4U 0142+61 also implies the presence of large toroidal 
magnetic fields and the possibility of prolate-shaped neutron 
star (Makishima et al. 2014). By considering possible growth 
of magnetic fields of magnetars during proto-magnetar phases, 
strong differential rotation within proto-magnetars would am¬ 
plify their toroidal magnetic fields (Duncan & Thompson 1992; 
Spruit 2009). Therefore, it would be natural that some magne¬ 
tars sustain large toroidal magnetic fields inside. 

The large toroidal fields are required from the stability anal- 
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yses of magnetic fields. Stability analyses have shown that stars 
with purely poloidal fields or purely toroidal fields are unsta¬ 
ble (Markey & Tayler 1973; Tayler 1973). Stable magnetized 
stars should have both poloidal and toroidal magnetic fields. 
Moreover, the toroidal magnetic field strengths of the stable 
magnetized stars have been considered to be comparable with 
those of poloidal components (Tayler 1980). However, we have 
not yet known the exact stability condition and stable magnetic 
field configurations, because it is too difficult to carry out sta¬ 
bility analyses of stars with both poloidal and toroidal magnetic 
fields. 

Nevertheless, stabilities of magnetic fields have been stud¬ 
ied by performing dynamical simulations. Braithwaite & Spruit 
(2004) showed that twisted-torus magnetic field structures are 
stable magnetic field configurations on dynamical timescale. 
Stabilities of purely toroidal magnetic field configurations or 
purely poloidal magnetic field configurations have been stud¬ 
ied in the Newtonian framework (Braithwaite 2006, 2007) and 
in the full general relativistic framework (Kiuchi et al. 2011; 
Lasky et al. 2011; Ciolfi et al. 2011; Ciolfi & Rezzolla 2012). 
Braithwaite (2009) and Duez et al. (2010) have found a stabil¬ 
ity criterion of the twisted-torus magnetic fields. It could be 
expressed as 


M Mp 
! \W\ < M 


< 0 . 8 , 


( 1 ) 


where M/\W\ is the ratio of the total magnetic energy to the 
gravitational energy. M p /M is the ratio of the poloidal mag¬ 
netic field energy to the total magnetic field energy, a is a cer¬ 
tain dimensionless factor of order of 10 for main-sequence stars 
and of order 10 3 for neutron stars. The ratio of M/\W\ is a 
small value (~ 10 -5 ) even for magnetars. Therefore, the crite¬ 
rion becomes 


0-2 < —tt 0.99, (2) 

M 

where Mt is the toroidal magnetic field energy. Stellar mag¬ 
netic fields would be stable even for toroidal magnetic field 
dominated configurations. Therefore, it is very natural that the 
toroidal magnetic field strength of the stable stationary magne¬ 
tized stars are comparable with or larger than those of poloidal 
component. 

Until recently, however, almost all numerically obtained 
equilibrium configurations for stationary and axisymmetric 
stars have only small fractions of toroidal magnetic fields, typ¬ 
ically Mt/M. ~ 0.01, even for twisted-torus magnetic field 
configurations in the Newtonian gravity (Tomimura & Eriguchi 
2005; Yoshida & Eriguchi 2006; Yoshida et al. 2006; Lander & 
Jones 2009; Lander et al. 2012; Fujisawa et al. 2012; Lander 
2013, 2014; Bera & Bhattacharya 2014; Armaza et al. 2014), 
in general relativistic perturbative solutions (Ciolfi et al. 2009; 
Ciolfi et al. 2010), and general relativistic non-perturbative so¬ 
lutions under both simplified relativistic gravity (Pili et al. 2014) 


and fully relativistic gravity (Uryu et al. 2014). All of them do 
not satisfy the stability criterion mentioned above. 

On the other hand, there appeared several works which 
have successfully obtained the stationary states with strong 
toroidal magnetic fields by applying special boundary condi¬ 
tions. Glampedakis et al. (2012) obtained strong toroidal mag¬ 
netic field models imposing surface currents on the stellar sur¬ 
face as their boundary condition. Duez & Mathis (2010) im¬ 
posed the boundary condition that the magnetic flux on the stel¬ 
lar surface should vanish. Since the magnetic fluxes of their 
models are zero on the stellar surfaces, all the magnetic field 
lines are confined within the stellar surfaces. They obtained 
configurations with strong toroidal magnetic fields which are 
essentially the same as those of classical works by Prendergast 
(1956), Woltjer (1959a, 1959b, 1960) and Wentzel (1960, 1961) 
and recent general relativistic works by Ioka & Sasaki (2004) 
and Yoshida et al. (2012). 

It is very recent that Fujisawa & Eriguchi (2013) have found 
and shown that the strong toroidal magnetic fields within the 
stars require the non force-free current or surface current which 
flows in the opposite direction with respect to the bulk current 
within the star. Such oppositely flowing currents can sustain 
large toroidal magnetic fields in magnetized stars. It is also very 
recent that Ciolfi & Rezzolla (2013) have obtained stationary 
states of twisted-torus magnetic field structures with very strong 
toroidal magnetic fields using a special choice for the toroidal 
current. Their toroidal currents contain oppositely flowing cur¬ 
rent components and result in the large toroidal magnetic fields, 
although their paper does not explain the physical meanings for 
appearances of such oppositely flowing toroidal currents. They 
also did not show clear conditions for the appearance of the 
toroidal magnetic field dominated stars. 

On the other hand, strong poloidal magnetic fields make stel¬ 
lar shape oblate one (e.g. Tomimura & Eriguchi 2005), but the 
strong toroidal magnetic field tends to stellar shape prolate one 
(Haskell et al. 2008; Kiuchi & Yoshida 2008; Lander & Jones 
2009; Ciolfi & Rezzolla 2013) Since the Lorentz force is a cross 
product of the current density and magnetic field, such Lorentz 
force requires a special relation between the magnetic fields and 
the current density. At the same time, the large toroidal mag¬ 
netic fields in stars also need a special relation between current 
density and Lorentz force. The oppositely flowing toroidal cur¬ 
rent density is a key to reveal these relations and a condition for 
appearance of the toroidal magnetic field dominated stars. 

We analyze magnetic field configurations and consider the 
special relations in this paper. We find a condition for the ap¬ 
pearance of the toroidal magnetic field dominated stars, which 
was not described in our previous work (Fujisawa & Eriguchi 
2013). In order to show the relations and condition clearly, sim¬ 
plified analytical models are solved and we show examples of 
the prolate configurations and the large toroidal magnetic fields 
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within stars. This paper is organized as follows. The formu¬ 
lation and basic equations are shown in Sec. 2. We present 
analytic solutions and the relations. We also explain the im¬ 
portant role of oppositely flowing components of the k currents 
for appearance of the prolate shapes and the presence of the 
large toroidal magnetic fields using the relations in Sec. 3. 
Discussion and conclusions follow in Sec. 4. In Appendix 1 
the deformation of stellar shape and the gravitational potential 
perturbation are briefly summarized. 


2 Formulation and basic equations 


Stationary and axisymmetric magnetized barotropic stars with¬ 
out rotation and meridional flows are analyzed in this paper. 

Some authors have claimed that there do not exist dynam¬ 
ically stable barotropic magnetized stars (e.g. Mitchell et al. 
2015 ), but in our opinion, their arguments should be ap¬ 
plied only to isentropic barotropes. The magnetized barotropic 
stars are well defined concepts apart from the thermal sta¬ 
bility or convective stability due to the entropy distributions 
and/or due to the chemical composition distributions. Thus we 
will investigate mechanical equilibrium states of traditionally 
defined barotropes (e.g., Chandrasekhar & Prendergast 1956; 
Prendergast 1956) in this paper. 

Since for such configurations the basic equations and basic 
relations are shown, e.g., in Tomimura & Eriguchi (2005) and 
Fujisawa & Eriguchi (2013), we show the basic equations and 
basic relations briefly. 

The stationary condition for the configurations mentioned 
above can be expressed as: 

J 7 = -^ + J l*{*)d* + C, (3) 


where p, p, <j> g and C are the density, the pressure, the gravita¬ 
tional potential of the star and an integral constant, respectively, 
ip is the magnetic flux function, p is an arbitrary function of VP. 
The magnetic flux is governed by 

A*lp = — 4 - 7 rrsin# —. (4) 

c 

where 

+ (5) 

ydr 2 r 2 d6 2 r 2 sin# dd J 

and j,p is a ^-component, i.e. the toroidal component, of the 
current density. The spherical coordinates ( r,6,ip ) are used. 

From the integrability condition of the equation of motion, 
the axisymmetry and the stationarity, the following relations are 
derived: 


i 

c 


1 dn 
47t d'P 


B + prsin#/r('P)e v 


( 6 ) 


« = *(*), (7) 

where j and B are the current density and the magnetic field, 


respectively, and k is another arbitrary function of \P. It would 
be helpful to note that the above relation for k was found by 
Mestel (1961) and Roxburgh (1966). Although re('P) is exactly 
a function of the magnetic flux function only in stationary and 
axisymmetric system (Braithwaite 2009), Braithwaite (2008) 
showed that the function k(\P) during the dynamical evolution 
of magnetized configurations is nearly conserved even for non- 
axisymmetric systems. 

Since the function k and the ^-component of the magnetic 
field is related as 

k(\P) = rsmdB v , ( 8 ) 


(9) 


the toroidal current density can be expressed as: 

j v 1 k('PW('P) . , . 

— = -- . n + prsm6p{V). 

c 47t rsm# 

Under our assumption that the magnetic field energy is small 
compared to the gravitational energy (M/\W\ < 10~ 5 ) in this 
paper, the influence of the magnetic fields can be treated as a 
small perturbation to a spherical star. Therefore, we assume 
that the stellar configurations are sphere and that the density 
profile depends only on r, i.e. p = p(r). For such situations, we 
can obtain analytical solutions easily. Noted that self-consistent 
approaches such as Tomimura & Eriguchi (2005), would re¬ 
veal some differences in the magnetic field solutions. In self- 
consistent approaches, we need to calculate both magnetic fields 
and matter equations iteratively. The stellar shape is no longer 
sphere and the stellar configuration affects the magnetic field 
configuration. However, our result in this paper is simple and 
might be important for both perturbative and self-consistent ap¬ 
proaches. 


3 Spherical models with weak magnetic 
fields 

Our aim in this paper is investigating the condition for appear¬ 
ance of the toroidal magnetic field dominated star analytically. 
Noted that our solutions of themselves are classical and not 
new ones, but we use the solutions in order to show the special 
condition clearly. 

3.1 Green’s function approach and analytic solutions 

We follow mostly the formulation of the classical works 
(Chandrasekhar & Prendergast 1956; Prendergast 1956; Woltjer 
1959a; Woltjer 1959b, 1960; Wentzel 1960, 1961) and the re¬ 
cent analytical works (Broderick & Narayan 2008; Duez & 
Mathis 2010; Fujisawa & Eriguchi 2013). In order to obtain 
analytical solutions, we choose the functional forms as follows: 

pW = Po, (10) 


k(^) = k 0 'I' , 


(ID 
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where po and ko are two constants. It should be noted that these 
functional forms always lead to non-zero surface currents unless 
magnetic fields are confined inside the star. The surface current 
induces a Lorentz force at the stellar surface (Lander & Jones 
2012). It would be unphysical because the Lorentz force need 
to be balanced by other physics such as a crust of the neutron 
star (e.g. Fujisawa & Kisaka 2014). The models relying on a 
surface current might not be physically realistic. We emphasize 
that we are not asserting that surface currents themselves are 
necessarily significant in real stars. The surface current simply 
provides a mathematically convenient way of describing analyt¬ 
ical solutions easily. 

By using these functional forms, the toroidal current density 
can be expressed as: 

= J- 'r, + Pop{r)rsind. (12) 

c 4-7T r sm 0 

We name the first term k current j Kv (force-free) term and the 
second term p current j£ (non force-free) term, respectively 
(Fujisawa & Eriguchi 2013). Then, the equation for the mag¬ 
netic flux becomes as follows: 

A*'!' + KqIF = — 4npop(r)r 2 sin 2 9. (13) 


It should be noted that this is a linear equation for the magnetic 
flux function with the inhomogeneous term which contains mul¬ 
tipoles less than the quadrupole. If we impose the boundary 
condition = 0 at the center of the star, 4/ is described as fol¬ 
lows (Duez & Mathis 2010; Fujisawa & Eriguchi 2013): 


T 


sin 2 9 


= Khiorji ( K 0 r ) 


- 47T/io«:oj rji ( K 0 r ) J y 1 (k 0 r') p(r')r' 3 dr 


+ ryi {n 0 r) / ji (ftor') p(r')r' 3 dr' >. 


(14) 


where we set the stellar radius r B = 1 in this paper, j i and yi are 
the spherical Bessel functions of the first kind and the second 
kind, respectively and K is a coefficient which is determined by 
a boundary condition of 4/ at the surface. According to the In¬ 
dependency of the inhomogeneous term, we search for solutions 
of the following form: 


a{r)sin8 = 'l>{r,e). 


(15) 


Therefore we obtain the solution for a(r ) by imposing the 
boundary condition at the surface and integrating equation (14). 

In this paper, we treat spherical polytropes with the poly¬ 
tropic indices N = 0 and N = 1. As for the configurations of 
the magnetic fields, we choose two types: (1) closed field mod¬ 
els (e.g. Duez & Mathis 2010) and (2) open field models (e.g. 
Broderick & Narayan 2008). For closed field models, since all 
magnetic field lines are closed and confined within the star, the 
magnetic flux must vanish at the stellar surface as follows: 


a(r s ) = 0 . 


(16) 


For open field models, since the poloidal magnetic field lines 
must continue smoothly through the stellar surfaces into the out¬ 
side. the boundary condition can be expressed as: 


a{r B ) 


da(r) 

dr 


The density profiles are 


(17) 


P(r) = Pc 


(18) 


for N = 0 polytrope and 

, , sin(7rr) 

P(r) = pc --—• 


(19) 


for N = 1 polytrope and p c is the central density. We can obtain 
four different analytical solutions according to four different sit¬ 
uations. We name them as a 0 c(r) (N = 0 with closed fields), 
ao o(r) (N = 0 with open fields), a\c{f) (N = 1 with closed 
fields) and aio(r) (N = 1 with open fields). 

Since the poloidal magnetic field lines are continuous 
smoothly at the surfaces for open field models, their external 
solutions a ex (r) must be expressed as: 


a ex (r) = _ ( 20) 

r 

It should be noted that the poloidal magnetic fields for closed 
field models are discontinuous at the surfaces except for solu¬ 
tions with special values of «o, i.e. eigen solutions with cor¬ 
responding eigenvalues (Broderick & Narayan 2008; Duez & 
Mathis 2010; Fujisawa & Eriguchi 2013). Therefore, these 
non-eigen configurations have toroidal surface currents. On the 
other hand, the toroidal magnetic fields for open field models 
are always discontinuous because of the choice of the functional 
form of k (equation 11). It implies that the open field configu¬ 
rations have non-zero poloidal surface currents. 

For non-eigen solutions, the toroidal surface current density 
can be expressed as follows: 


jip,sur(P) _ 1 

47T 




1 

( 

av? ir *\ 

47rr s sin 6 

^ dr 

dr ) 


sin# 

( da ex 

da in \ 

47rr s 

K dr 

dr J 


= jo sin d, 


( 21 ) 


where superscript 171 denotes an internal solution and jo is a 
coefficient of the surface current density. 

Analytic solutions are obtained by fixing a boundary condi¬ 
tion and integrating equation (14). Four different inner solutions 
(0<r<l) can be obtained according to four different situations 
as follows: 
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k 0 = 4.0 



*0 = 7-0 



Fig- 1- Distributions of the toroidal current density normalized by the maximum strength of I'l'maxl are shown along the equatorial plane. Curves with different 
types denote the behaviors of the total toroidal current density, 3 v /c, (thick solid line), the toroidal k 0 current density (thin solid line) and the toroidal mo current 
density (thin dotted line). We set mo = -1 in order to plot these distributions. Left panels show the profiles of solution a 10 with k 0 = 1.0 and 4.0 and right 
panels show those of solution aic with k 0 = 2.0 and 7.0. 


aoc(r) = 4npoPc 


sin(Kor) — Korcos(Kor) 
rtt q (sin k o — k q cos kq ) 



( 22 ) 


at «o = 7r (Haskell et al. 2008), while the solution aoc is not 
singular at ko = n (Fujisawa & Eriguchi 2013). 


do o(r) = 4-rrpoPc 


3{sin(/tof) — Korcos(Kor)} r 2 
tKq sin k o Kg 


(23) 


aic(r) = 


Pope 


87r{sin(Kor) — Kor cos(Kor)} 
sin kq — Kq cos Kq 


- { (4kq — 47T 2 ) r 2 + 8} sin(7rr) + 8irr c.os(jrr) 


, (24) 


aio(r) 


Pope 

r(* 2 o-i r 2 ) 2 
- {(4«o - 


(47t 3 — 47TK§){sin(Kor) — Kor cos(Kor)} 
Kg sin Ko 

47r 2 )r 2 +8} sin(7rr) + 8nrcos(irr) .(25) 


The open field models (aoo and aio) continue to the external 
solutions (r > 1) expressed by equation (20). Here it would be 
helpful to explain several different kinds of characteristic solu¬ 
tions. 

First, for aic(r) solutions there appears a singular solution 


Second, although most solutions are accompanied by surface 
currents, some special solutions have no surface currents. We 
call such solutions without surface currents as eigen solutions 
and the values of ko as eigenvalues. 

Third, there appear many eigen solutions as the value of 
k o exceeds the first eigenvalue. We call those eigen solu¬ 
tions as higher-order eigen solutions (see figures in Broderick 
& Narayan 2008; Duez & Mathis 2010; Yoshida et al. 2012). 
Those solutions appear when the value of ko exceeds the first 
eigenvalue of ko for each situation. 

Fourth, special solutions with different polytropic indices 
come to coincide with each other. In other words, those so¬ 
lutions do not depend on the matter distributions. As seen from 
the expression for the current density, the contribution from the 
p current term needs to disappear. It implies that those solutions 
are determined only by the k current. Since the k currents do 
not contribute to the Lorentz force, these solutions are called as 
the force-free solutions (Wentzel 1961). The force-free solution 
is expressed by the following form: 
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a //( r ) = -{sin(«oJ") — Korcos(Kor)}. (26) 

«o r 

The solution becomes force-free when ko ~ 4.49 and 7.73 for 
closed models and when k 0 = 7r and 2-7T for open field models. 
(Broderick & Narayan 2008; Fujisawa & Eriguchi 2013). 

The toroidal surface current in equation (21) vanishes when 
the ko is eigenvalue, i.e., for eigen solutions. The lowest eigen¬ 
values are k 0 ~ 5.76 for a 0 o, ~ 7.42 for ai c , ~ 5.76 for a 0 o 
and ko ~ 4.66 for aio. Hereafter, we focus on solutions with 
ko less than the lowest eigenvalue. However, our analyses and 
results could be general and would be valid even when the con¬ 
figurations are higher-order eigen solutions. 

In figure 1, distributions of the normalized j v /c (thick solid 
line), the k current (thin solid line) and the p current term 
(dashed line) along the equatorial plane are shown for solutions 
of aio (with ko = 1.0,4.0, smaller and larger than force-free 
ko, respectively) and me (with ko = 2.0,7.0, smaller and larger 
than force-free ko, respectively.) We have fixed po = — 1 fol¬ 
lowing Fujisawa & Eriguchi (2013) in order to plot these curves. 

As seen in upper panels in figure 1, directions (signs) of the 
p current, i.e. non force-free current, and the k current, i.e. 
force-free current, are the same for solutions with smaller ko- 
By contrast, for solutions with larger kq (lower panels in figure 
1), the p current flows oppositely to the k current (Fujisawa & 
Eriguchi 2013). Moreover, most of the j^/c (thick solid line) 
for solutions with ko = 4.0 and k o = 7.0 flows oppositely against 
the corresponding p current. Since the sign of the total toroidal 
current determines the sign of the magnetic flux function (see 
equation 4), this implies that the sign of po 4/ for the whole inte¬ 
rior region changes from po^ > 0 to po^fr < 0 at the force-free 
solutions. We calculate many solutions and confirm that the the 
sign of po'i for the whole interior region changes at the force- 
free solution. We call the current distribution for which po^ < 0 
oppositely flowing current. 

On the other hand, the surface toroidal currents in the closed 
field models are always oppositely flowing to the total toroidal 
currents because of the zero-flux boundary condition equation 
(16) and the form of the surface current equation (21). 


cally in terms of the current density. This is the reason why we 
consider both magnetic fields and current density in this paper. 
As we have seen in the Sec.3.1, the oppositely flowing toroidal 
current density ( po^ < 0 ) plays a key role for appearance of 
the large toroidal magnetic fields. The direction of the toroidal 
current seems to relate to the stellar deformations because the 
Lorentz force is a cross product of current density and magnetic 
held. We consider the relation between the toroidal current and 
the poloidal deformation of stars in this subsection. 

In our analytic models, the Lorentz force L is expressed us¬ 
ing the arbitrary function p('V) as 


L = 


(^xsj = pV J p{\ V) d® = p/r(®)V® 

da . 2 n „ a . 

= ppo — sm 9e r + 2ppo — smftcost/eg. 
dr r 


(27) 


Following Haskell et al. (2008), we consider the stellar 
quadrupole deformations of N = 1 polytropic stars. Haskell 
et al. (2008) calculated magnetic deformations of polytorpic 
magnetized star with poloidal and toroidal magnetic fields. 
Although they derived the general forms of the deformations 
(equation 64, 65 & 67 in their paper), they did not show the an¬ 
alytical expressions of them explicitly. They displayed only a 
few numerical results in Tab. 1 in their paper. By contrast, we 
show the analytical solutions of the deformation in order to in¬ 
vestigate the condition for appearance of the toroidal magnetic 
field dominated star. 

We assume that the influence of the magnetic fields to the 
stellar structures are small and that their effects can be treated 
perturbatively. Due to the effects of the magnetic fields, a cer¬ 
tain physical quantity X (r, 9) is assumed to be expressed as 


X{r,9) = X{r) + ^5X {n \r)P n (cos 9 ), (28) 

n=0 

where SX^ denotes a small change of order 0(B 2 ) of the 
quantity X due to the Lorentz force. The angular dependencies 
are treated by the Legendre polynomial expansions and the co¬ 
efficient of each Legendre polynomial is expressed as SX ^ (r). 
This expansion is also applied to the Lorentz force as follows: 


3.2 Deep relation between the toroidal current and 
the poloidal deformations of stars 

As the many previous works pointed out, the toroidal magnetic 
fields tend to deform stellar shapes prolate, while the poloidal 
magnetic fields tend to deform them oblate (Wentzel 1960, 
1961;Ostriker & Gunn 1969; Mestel & Takhar 1972). These 
studies used only the magnetic fields in their formulations. The 
ideal MHD system can be described by using only magnetic 
fields and one does not need to mention the electrical current 
density at all. In contrast, we consider both magnetic fields 
and current density in our calculation. Although these two ap¬ 
proaches are equivalent, it is easier to interpret results physi¬ 


L(r,6O = ^L (n) (r)P„(cos0). (29) 

71=0 


From the perturbed equilibrium condition equations, the follow¬ 
ing relations can be derived: 


dSp <"> , JStfr^ , r „ (n) # s r(B) 

-;-M- - - \-0p —— - , 

dr dr dr 


(30) 


Sp' 71 ' 1 + pScjtp = rL. 


(31) 


Since we are interested in the quadruple deformation, we con¬ 
sider only n = 2 components of Lorentz force as follows: 

( 2 ) _ 2pp 0 da(r ) 

±J r —- 


3 dr 
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( 2 ) _ 2,pp 0 ffl(r) 


rW _ 

L e — 


l ( 2) = ^(2) _ rf(^4 2) ) 

r dr 

2po dp 

-a(r). 


(32) 


3 dr 

The change of the stellar surface to the order of the quadrupole 
term can be expressed as 

rd(0 ) = r s {l + eP 2 (cos#)} = r s jl + ^(3 cos 2 # — 1) j , (33) 

where ra(9) denotes the deformed surface radius and e is a 
small quantity which represents the fraction of the stellar sur¬ 
face along the pole. Following this expression, the stellar shape 
is prolate for e > 0 and oblate for s < 0. 


3.2.1 Deformation of N ^ 0 polytrope 

Using these equations, the quadrupole change of the density is 
described by: 

„(2) _ ( dp x A 2) , r(2)\ ( d <t>9 N _1 




(34) 


Since the surface of the deformed star is defined by a set of 
points where the pressure vanishes, i.e. 


p{r d {9)) = 8p(r s ) +er s P 2 (cos#) 

ar 


= 0 


we can derive 


_ dpo 
dr 


+ 8p 


( 2 ) 


= 0 


(35) 


(36) 


for polytropes with N ^ 0. For N = 0 polytrope, this equation is 
reduced to trivial relation 0 = 0 and so we will treat N = 0 poly¬ 
trope differently as will be shown in the next section. Therefore, 
the quadrupole surface deformation e for IV ^ 0 is obtained by 

' dp\~ 1 5pW 


( 2 )' 


(37) 


It is clearly seen that, since < o at the surface, the stellar 
deformation is prolate for 8p^ > 0 and oblate for 8p^ < 0. In 
our situation, the explicit form of <5p (2) can be expressed as 


and e for N ^ 0 polytropes becomes as 

_ / 2pp 


(38) 


e = - ^-^a(r s ) + 8^ (r)) 


#g 2) 

dr 




(39) 


As shown in Appendix 1 the gravitational change for N = 1 
polytrope can be obtained as 




5^(x) 


dx 


h(x). 


(40) 


^g 2) O) = 


fW(tt) 


3 dF (p) 

7T 4 dx 


(41) 


where = 3/n 2 is used. Here the function F (:P> (x) is de- 
hned in Appendix 1. This is a analytic solution of the deforma¬ 
tion. 

Since the expression for the function F is so compli¬ 
cated, it is not clearly seen the sign of the quantity (8rj£p (r s ) + 
2poa{r s )/3) which determines the sign of the quantity e. In hg- 
ure 2 we show the behaviors of —8rj^g\r B ) and —2poa(r s )/3 
against the value of kq. As we have seen, the sign of po'b 
changes at the force-free solution ko ~ 4.49 for closed model 
and ~ 7r for open model. As seen in this figure, the shape change 
from the effect due to the gravitational change is the same as that 
from the Lorentz term. Thus the sign of the quantity £ is essen¬ 
tially determined by the sign of the Lorentz term, i.e., the sign 
of the quantity poa(r B ). Since p(r)(dp/dr)~ 1 < 0, the stellar 
shape is oblate for /to 4'(r, #) > 0 for the whole interior region 
and prolate for po 4' (r, 9) < 0 for the whole interior region as far 
as the global poloidal magnetic field is dipole. Therefore, the di¬ 
rection of the deformation by Lorentz force is determined by the 
direction of the non-force free current p current (equation 12). 
If the ^current flows oppositely to the magnetic flux (po’L < 0), 
the stellar shape is prolate. If the p current flows same direc¬ 
tion (po'L > 0), the stellar shape becomes oblate one. This is a 
deep relation between the direction of the toroidal current and 
the poloidal deformations of stars. 

In figure 3, the contours of 4> (dashed curves) and the direc¬ 
tions of Lorentz force vectors (arrows) are displayed. It should 
be noted that directions of the Lorentz forces are totally oppo¬ 
site between models with /.to4* (r B ,9) > 0 (ko = 1.0 and 2.0) and 
those with po^rs,#) < 0 (kq = 4.0 and 7.0) 


3.2.2 Deformation for N = 0 polytrope 

For N = 0 polytrope, the gravitational change and the shape 
change are written as follows as shown in Appendix 2: 


S0g 2) = --7T Gpoer 2 , 


rg 

and 


{-^Gpor 2 - (-^TrGpor 2 ) j ^ a{r s ) 


5p 0 


■a(r s ). 


AnGpo 

Here we use the stationary condition 

8p^ + p 0 8^ =rL^ , 


(42) 


(43) 


(44) 


and the surface condition equation (35). 

Thus the sign of the quantity e is exactly determined by the 
sign of the Lorentz term, i.e., the sign of the quantity poa(r s ). 
The stellar shape is oblate for /to<E'(r,#) > 0 for the whole inte¬ 
rior region and prolate for po4/(r,#) < 0 for the whole interior 
region as far as the global poloidal magnetic field is dipole. The 


Thus for x = 7r, i.e. on the surface, 
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Fig. 2. The values of - 2 p 0 a(x = 7t)/3 (thin solid line) and — (x = -n-) (thin dashed line) in closed field model (left panel) and open field model (right 
panel) are plotted. The thick vertical lines denotes force-free limit. The toroidal current densities consist of oppositely flowing flows beyond the dashed thick 
vertical lines. We set yi 0 = — l and p c = l in order to plot these graphs. 


Kg-1.0 


Kq-4.0 


k 0 =2.0 


k 0 =7.0 



R 




Fig. 3. Poloidal magnetic field structures (dashed curves) and Lorentz force vector fields (arrows) for the open field models (k 0 = 1.0, /to = 4.0) and the 
closed field models (k 0 = 2.0 and k 0 = 7.0) are displayed.Vectors only show their directions but are not scaled to their absolute values. 


relation between the direction of the p current and the poloidal 
deformations of stars is still valid in this case. 

3.3 Deep relation between the toroidal current and 
the strong toroidal magnetic fields 

We have found a relation between the oppositely flowing 
toroidal current density and the Lorentz force in the previous 
subsection. We consider a relation between the oppositely flow¬ 
ing toroidal current and the strong toroidal magnetic fields in 
this subsection. 

In figure 4, the ratio of the toroidal magnetic field energy 
Mt to the total magnetic field energy M — M v + Mt of each 
model is plotted for different situations. The solution becomes 
force-free at the point denoted by the vertical solid lines. The 


dashed vertical lines denote the critical values beyond which 
there arise oppositely flowing k current structures dominated. 

As seen in figure 4, N = 0 solutions and N = 1 solutions 
cross at ko ~ 4.49 and 7.73 for closed field models and ko = n 
and 27 t for open field models, because the solutions at these 
points are force-free solutions as mentioned before. The energy 
ratio is Mt/M ~ 0.5 when the solutions are the first force- 
free configurations. Therefore the solutions are divided into two 
types at the force-free solution. The solution whose /to value is 
smaller than force-free /to is poloidal dominant configuration, 
while the solution with larger /to is toroidal dominant config¬ 
uration. Since the sign of /to 4' changes at the force-free so¬ 
lution, the solution is poloidal dominant for po'hir, 6) < 0 for 
the whole interior region (oppositely flowing current) and the 
solution is toroidal dominant for /to 4 '(r,6) > 0 for the whole 
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Closed field model Open field model 




K o 


K o 


Fig. 4. Energy ratio M t /M is plotted against the value of k 0 . Closed (left panel) and open (right panel) field solutions are shown. The solid and dashed 
curves denote N = l and N = 0 solutions, respectively. The vertical solid lines mean force-free solutions: Closed force-free solutions appear at k 0 ~ 4.49 and 
k 0 ~ 7.73 and open force-free solutions at k 0 = ir and k 0 = lix. The toroidal current densities are composed of two oppositely flowing components beyond 
the vertical dashed lines: k 0 ~ 5.76 for the a 0c solution (dashed curve in left panel), k 0 ~ 7.42 for the a lc solution (solid curve in left panel), k 0 ~ 5.76 for 
the a 0o solution (dashed curve in right panel) and k 0 ~ 4.66 for the a l0 (solid curve in right panel). The open circle in the left panel denotes the singular 
solution for a lc (r). 


interior region. The oppositely flowing non-force free current 
(/jo it (r, 9) < 0 for the whole interior region) is required for large 
toroidal magnetic fields. This is a relation between the toroidal 
current density and the toroidal magnetic field. 

3.4 A situation for appearance of toroidal magnetic 
field dominated configurations 

As we have shown in previous parts in this paper, there are two 
deep relations between toroidal current, poloidal deformation 
and strong toroidal magnetic field. One is a relation between 
the toroidal current and the poloidal deformation of stars in Sec. 
3.2. The other is a relation between the toroidal current and the 
strong toroidal magnetic fields. The important finding in this pa¬ 
per is that the appearance of oppositely flowing non force-free 
current which fulfills the condition (/io'D < 0) changes the stel¬ 
lar shape prolate and makes the toroidal magnetic fields toroidal 
dominant. Therefore, a well-known relation between toroidal 
dominant magnetic fields and prolate shapes requires the oppo¬ 
sitely flowing non force-free toroidal current density. Although 
our result is very simple and natural, nobody have described ex¬ 
plicitly that the oppositely flowing non force-free current den¬ 
sity makes the stellar shape prolate. It might be because almost 
all previous studies treated only magnetic fields and did not pay 
special attention to current density. 

Consequently, we can conclude that a condition for appear¬ 
ance of prolate configurations and the toroidal magnetic field 
dominated configurations is that the arbitrary function sat¬ 
isfies the following condition: 

p(^)d\\f < 0 , (45) 

for the whole interior region when the functional forms are 


equations (10) and (11). Although, exactly speaking, these anal¬ 
yses and conditions are valid within the present parameter set¬ 
tings, our results would be useful for more general situations. 
This might be naively seen from the contribution of the term 
f pd^> in the stationary condition (equation 3). If this term is 
negative, it implies that the action of the Lorentz term is oppo¬ 
site to that of the centrifugal force which is expressed by the 
term f fl(R) 2 RdR and is always positive. In other words, the 
magnetic forces or Lorentz forces act as if they are the ’anti’- 
centrifugal forces and therefore shapes of stationary configu¬ 
rations become prolate (see also calculations in Fujisawa & 
Eriguchi 2014). 

Although the condition presented in this paper might not 
be always correct, we could obtain the large toroidal magnetic 
fields by employing this criterion for more complicated calcu¬ 
lations. 

4 Discussion and Summary 

4.1 Physical reason for the necessity of appearance 
of k currents to realize prolate configurations 

In order to get configurations with prolate shapes, we need to 
include the ’anti’-centrifugal effects or ’anti’-centrifugal poten¬ 
tials. As is easily understood, the anti-centrifugal potentials 
should behave as decreasing functions from the symmetric axis, 
or, at least, they must contain decreasing branches which cover 
wide enough regions to result in effectively anti-centrifugal ac¬ 
tions. 

For our formulation, the following property is commonly 
found: 

p > 0 ->• j£> 0 


d/ > 0 -> / pMt > 0 , (46) 
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Fig. 5. Distributions of the Mr (dashed line) and d J fj,d^/dR (solid line) of closed field solutions are plotted. The left panel shows the distributions with 
k 0 = 2.0 and /xo = l.O and the right panel shows those with k 0 = 7.0 and /j, 0 = —1.0. 


and 

P < 0 -»■ j£ < 0 


® < 0 


pdfS > 0 . (47) 


In addition to these behaviors, for > 0 configurations, the 
magnetic flux functions increase to the maximum values as the 
distance from the axis and turn to decrease beyond the maxi¬ 
mum point as follows: 


d f pdS 
d R 

d J pdS 

d R 


>0 

<0 


for R< R max 5 
for R> R 

max 5 


(48) 

(49) 


where R max is the location of the maximum point of the mag¬ 
netic flux function 4' (left panel in figure 5). 

For < 0 configurations, the magnetic flux functions de¬ 
crease to the minimum values as the distance from the axis and 
turn to increase beyond the minimum point: 


d J pd<S 
d R 

d J pdS 
dR 


>0 

<0 


for R fini i N , 
for R > Rmin , 


(50) 

(51) 


where R m i n is the location of the minimum point of the mag¬ 
netic flux function 4'. 

Although there exist decreasing branches for both situations, 
these decreasing branches cannot overcome the centrifugal ef¬ 
fects due to the increasing branches. Therefore, the global con¬ 
figurations with purely ip-currents would become oblate shapes. 

From this consideration, the anti-centrifugal forces could be 
realized if the following (necessary) condition is fulfilled: 


p > 0 AND j v < 0 AND ^ < 0 AND 


pd'S < 0 ,(52) 


or 


p< 0 AND j v > 0 AND ® > 0 AND / pdS < 0 .(53) 


These conditions could be realized only by including the n- 
currents so that the following conditions are satisfied: 


> 0 , j v =jv+j% < 0 AND $ < 0 , (54) 

for R < R min , (55) 

for R > -R m in , (56) 

3v < 0 , 3v = &+& > 0 AND S > 0 , (57) 

for R < Rmax , (58) 

for R > R max . (59) 

in figure 5 shows the distributions of 'll and 

= 7.0 and p = — 1. As seen in figures 1 and 
5, the conditions of the above-mentioned are satisfied undoubt¬ 
edly. Therefore, the appearance of k currents j" which are op¬ 
positely flowing with respect to the p currents and at the 
same time whose magnitudes are large enough to overcome the 
p currents are required to realize prolate shapes. 

4.2 Twisted-torus configurations with large toroidal 
magnetic fields 

Almost all previously carried out investigations for magnetized 
equilibrium states having twisted-torus magnetic fields had 
failed to obtain toroidal magnetic field dominated (Mt > M p ) 
models. We have found that most models of their works do not 
satisfy the condition of equation (45) and the magnetized stellar 
shapes are oblate due to the p current term. The k term in those 
works has been chosen as follows: 

«(*) = K0(V& - ^ m ax) fcl+1 0(^ - 4) max) , (60) 

where, k\ is a constant and 0 is the Heaviside step function 
and T-max is the maximum value of on the last closed field 
line within the star. Since the current density of this functional 


P > 0 , j% < 0 . 

d J pd<S 
OR 

d f pdS 


dR 


<0 

>0 


p < 0 , j% > 0 . 

d f pd'S 
dR 

d f pd'S 


<0 


dR 


>0 
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form vanishes at the stellar surface, there exist no surface cur¬ 
rent and exterior current density. This functional form was used 
by Tomimura & Eriguchi (2005) for the first time and results 
in the twisted-torus configurations. The same choice for the k 
has been employed by many authors (e.g. Yoshida & Eriguchi 
2006; Yoshida et al. 2006; Kiuchi & Kotake 2008; Lander & 
Jones 2009; Ciolfi et al. 2009; Ciolfi et al. 2011; Fujisawa et al. 
2012;Glampedakis et al. 2012; Lander et al. 2012; Fujisawa 
et al. 2013; Fujisawa & Eriguchi 2013; Lander 2013, 2014). 
While the functional form p(^) = po (constant) has been used 
in many investigations, Fujisawa et al. (2012) and Fujisawa 
et al. (2013) used a different functional form as 

M (^)=/ro(4> + e) m , (61) 

where m and t are positive constants. They have obtained 
highly localized poloidal magnetic field configurations using 
this type of functional from. However, their works did not sat¬ 
isfy the condition of equation (45) and did not obtain models 
with large toroidal magnetic fields. 

Recently, Ciolfi & Rezzolla (2013) have adopted a pertur¬ 
bative approach and succeeded in obtaining magnetized equi¬ 
librium states with twisted-torus magnetic fields whose toroidal 
fields are large. Their functional form of k is 

«(*) = K 0 'I'(|>I'/^max| - 1)6 ( | V&/^ m ax | - 1). (62) 

On the other hand, the functional form of p is 

M®) = Co[(l - |^/^ ma x|) 2 0(l - | 4/ / 4/ max |) - k] 

+ Xo (63) 

where co, k(> 0) and Xo are constants. The toroidal magnetic 
field is confined within the last closed field line in these func¬ 
tional forms. Outside the toroidal magnetic field region, the 
function k vanishes and p becomes 

/r(®)= Co [(l-|^/^ m ax|) 2 -fc] ■ (64) 

Since the first term and the second term are positive and neg¬ 
ative, respectively, this function with larger k tends to satisfy 
the condition of equation (45). As they noted, larger values of 
k result in larger energy ratios Mt/M. As the value of k in¬ 
creases, the energy ratio Mt/M. increases and the stellar shape 
becomes more prolate in general (see Tab. 1 in Ciolfi & Rezzolla 
2013). However, they assumed that the magnetic field config¬ 
uration is purely dipole but their functional forms and toroidal 
current density distribution are far from dipole one (see bottom 
panels of figure 2 in Ciolfi & Rezzolla 2013). Non-perturbative 
studies with higher order components were unable to repro¬ 
duce their results and found contradictory results (Bucciantini 
et al. 2015). We need to calculate magnetic field configurations 
with higher order components for large toroidal models by us¬ 
ing non-perturbative methods in the future. 

The condition of equation (45) itself is valid when a star is 


barotropic. However, the relation between oppositely flowing 
toroidal current density and prolate share is very simple and 
natural when a star is non-barotropic. Therefore, this condi¬ 
tion is also useful for recent perturbative non-barotropic solu¬ 
tions (Mastrano et al. 2011; Mastrano & Melatos 2012; Akgiin 
et al. 2013; Yoshida 2013). We also need to investigate non- 
perturbative non-barotropic magnetized equilibrium states in 
the future. 

4.3 Summary 

In this paper, we have obtained four analytic solutions with both 
open and closed magnetic fields for spherical polytropes with 
weak magnetic fields. 

Using the obtained solutions we have discussed the situa¬ 
tions for which the prolate equilibrium states and the toroidal 
magnetic field dominated configurations appear. The main find¬ 
ing in this paper is that the appearance of the prolate shapes and 
the toroidal magnetic field dominated states are accompanied 
by the appearance of oppositely flowing k currents with respect 
to the p current. This situation seems to be related to the con¬ 
dition for the non force-free toroidal current contribution, i.e. 
f /r(4t)d4/, in the stationary state condition equation (3). 

Although the appearance of prolate shapes and the occur¬ 
rence of toroidal magnetic field dominated states cannot be de¬ 
fined quantitatively, the rough qualitative idea about them can 
be determined by checking the sign of the magnetic field poten¬ 
tial, i.e. the quantity f /r(4t)d4>. 

Of course, the analytic solutions obtained in this paper have 
been derived under very restricted assumptions. However, as 
explained in the Discussion, the concept of the ’anti’-centrifugal 
actions due to the magnetic potentials would be applied to more 
general situations for the magnetic fields. 
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Appendix 1 Change of the gravitational 
potential for N = l polytrope 

The gravitational potential perturbation for N = 1 polytrope is 
governed by the quadrupole component of Poisson’s equation 
under two boundary conditions is regular at r = 0 and 

continues the external solution smoothly at r = r 3 ): 

d 2 S(p (2) 


+ 2 d ^__6 ( 2)=4 (a) 

r dr r z 


(Al) 


dr 2 r dr r 2 

Considering the density perturbation expressed by equation 
(34), this equation can be rewritten as 


d 2 ScP ( g ) 
dr 2 


2 dScj) (2) 

+ -p- + 

r dr 


4 2) 


(tv 2 - 4-) ^9 


(A2) 


By introducing the new variable x = nr, the left-hand side of 
the equation is reduced to 

A( 2 ) 


dfs^g’ 2 d5(f)'. 


dx 2 


x dx 


_ 4nG / d<f> g \ 
7t 2 V dr ) 


( 2 ) 

9 


l/ 2 -* (r). 


(A3) 


The solution to this equation can be obtained by taking the 
boundary conditions into account as follows: 


( 2 )^ F^(x) 1 dF^(n) 


5^g J, (x) = 
where 


dx 




(A4) 


F (p \x) 

= - IftoAi ^ {67T 2 K 2 + (7T 2 -3 kq) kq* 2 } sin (—x) 
Kq(7H — Kq) 2 L l v ’ \n J 


3 

Kq 


- g/roA2 


( 2 

2\ 

2 2 ~l 

f /to Y 

(tt 

-« oj 

KqX j COS 

vV x )\ 


1 4 . . 1 

2 X SmX+ 6 




(M 

X 5 COS X 


(A5) 


Here the coefficients Ai and A 2 are defined as 


Ai = 


A 2 = 


8n 2 popc 


(kq — 7r 2 ) 2 (sin/to — «o cos/to) ’ 

4np 0 Pc 


(K 2 -n 2 ) 2 ’ 

for N = 1 closed configurations and 


Ai = — 


4n 2 pop c 


Kq(k,q — 7r 2 ) sin/to 


(A6) 


(A7) 


(A8) 
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An p 0 p c 


(A9) 


for N = 1 open configurations. 

Appendix 2 Surface change for N = 0 
polytrope 

The change of the gravitational potential due to the change of 
the surface, i.e. er s P 2 (cosd), can be obtained by 






